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Abstract: The possibility to compute nucleon electromagnetic form factors in the time-like region by analytic continuation 
of their space-like expressions has been explored in the framework of the Skyrme model. We have developed a 
procedure to solve analytically Fourier transforms of the nucleon electromagnetic current and hence to obtain 
form factors defined in all kinematical regions and fulfilling the first-principles requirements. The results are 
discussed and compared to data, both in space-like and time-like region. 
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1. Introduction 

Nucleon electromagnetic form factors [1] (FFs) are Lorentz scalar functions of the squared four-momentum trans¬ 
fer of the photon, q^, that parametrize those degrees of freedom of the nucleon electromagnetic current, which 
are not constrained by Lorentz and gauge invariance. 

They represent a unique source of information about the internal structure of nucleons. In particular, in the non- 
relativistic limit (low q^), FFs can be interpreted as the Fourier transforms of the electric charge and magnetic 
momentum spatial distributions of the nucleon. 

From the point of view of quantum held theory, being related to the electromagnetic current and hence only to 
the Born amplitude (one-photon exchange), FFs embody the resummation of all high order processes with two 
nucleons and one photon as external particles. 

Such high-order processes represent the connection with quantum chromodynamics (QCD) of FFs, that indeed 
could be described in terms of hadronic loops, involving virtual mesons and baryons. Due to the large number 
of hadron ’’species” to be accounted for and also to the unknown couplings among them, a direct calculations of 
FFs in the framework of QCD, especially in the low-q^ regime, is a very hard task. 

Nevertheless, interesting results have been obtained by lattice calculations [2] and effective model approaches, such 
as: chiral perturbation theory [3], chiral soliton models [4], large-Vc approximation [5] and holographic QCD [6]. 
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However, in the majority of these cases, the obtained FF descriptions are restricted to the only space-like (SL) 
region. 

In general, two kinds of FFs could be identified: SL FFs (SLFFs), related to the elastic scattering process 
eN —>■ eN {N stands for nucleon), which occurs with cf < Q (see, for instance, Refs. [7] and [8]), and time-like 
(TL) FFs (TLFFs), related to the annihilation processes e^e“ NN, where > {2Mn)^, Mn is the nucleon 
mass (for a review see Ref. [9] and references therein). 

The scattering and annihilation processes are related by crossing symmetry, which, considering only the Born ap¬ 
proximation, see fig. 1, implies that SLFFs and TLFFs represent values, for negative and positive respectively, 
of a unique function of , simply named FF. 

As a consequence, in order to understand the meaning of FFs, especially in the TL region, where the interpreta¬ 
tion in terms of Fourier transforms of spatial distributions fails, we must adopt descriptions or parametrizations 
defined in the whole kinematic region. Moreover, new FF data, coming from different experiments^, should help 
in shedding light especially in the more puzzling TL region. 



Figure 1. Feynman diagram in Born approximation for e^e NN and eN —>■ eN. 

symbolizes FFs. 


The solid disc at the nucleon vertex 


Various techniques and procedures have been proposed to develop such a SL-TL unified description of nucleon 
FFs. Many of them make use of dispersion relations, see for example Refs. [14-16], others suggest new models 
(for instance, in Ref. [17], a semi-phenomenological microscopic model is proposed) and some others use analytic 
continuation methods to extend, to all values of , parametrizations usually dehned only in the SL or TL re¬ 
gion [18]. 

We will expound here a procedure, originally formulated in Ref. [20], that allows to make the analytic contin¬ 
uation to the whole q^ complex plane of a parametrization of FFs, initially conceived for a particular reference 
frame. More in detail, nucleon FFs are computed in the Breit frame as Fourier transforms of the time and space 
components of the electromagnetic current. This particular representation is defined only in the SL region, i.e., 
the Fourier integrals, which depend on q^, converge only for SL four-momenta. However, if such a representation 
can be analytically solved, that is, the Fourier transforms are obtained as analytic functions of q^, instead of 
discrete numeric values at each four-momentum transfer, then the FF parametrizations should be valid in all 
non-singular points of the q^ complex plane. 

In particular, the nucleon electromagnetic current has been computed in the framework of the Skyrme 
model [19, 21], by solving numerically a set of non linear differential equations. The most relevant aspect of 
the procedure outlined here, consists in assigning to these numerical solutions, opportune analytic expressions, so 
that their Fourier transforms embody the properties required for FFs by analyticity and unitarity. 

The structure of the article is the following: in the second section we briefly introduce FFs in SL and TL regions 
and describe their analytic properties. In the third section we review the Skyrme model and calculate SLFFs. 
In the fourth section we illustrate the method of analytic continuation and the obtained results. In closing, we 
discuss the main issues of these results, also in comparison with all available FF data. Numerical computations 
have been performed by means of the software package Mathematica 7.0 [22]. 


^ BESIII [10] at BEPCII in Beijing, China; SND [11] and CMD3 [12] at VEPP-2M in Novosibirsk, Russia; PANDA at 
FAIR in Darmstadt, Germany [13]. 
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1.1. Space-like form factors 

In the scattering channel, fig. 1 vertical direction, the Feynman amplitude of the nucleon vertex, N ^ 'y*N, is 
parametrized as [23] 

(iV'(p')b'"(0)|iV(p)) = m(p') (9^)) “(P) ’ (1) 

where the four-momenta follow the labelling of fig. 1 and, {q^) and F 2 {q^) are the so-called Dirac and Pauli 
FFs (the ’’tilde” indicates their SL definition). They are Lorenz scalar functions and, as a consequence of the 
hermiticity of the current operator and the time reversal symmetry, are real for q^ <0. At q^ — 0 the Dirac 
FF is normalized to the nucleon charge Qiv, in units of the positron charge, while the Pauli FF is normalized to 
the anomalous magnetic moment kw, in units of the Bohr magneton, i.e., 

^i'^( 0 ) = Qiv, Fi^( 0 )=«iv. ( 2 ) 


In the special frame, called Breit frame, where there is no energy exchange, hence: p = {E, —ql2), p' = {E,q/2) 
and q = ( 0 ,^, the time and space components of the current expectation value, eq. ( 1 ), reduce to 


{N'{p')\J°{0)\N{p))^Fl'{q^) + ^Fi^iq^), 

{N'{p')\J{0)\N{p))=uip')yu{p)[Fr{q^) + Fi’{q^)) . 


(3) 


These combinations of the Dirac and Pauli FFs, representing the Fourier transforms of charge and magnetization 
spatial distributions of the nucleon, define the electric and magnetic Sachs FFs [24] 




(4) 


that, following eq. ( 2 ), are normalized at = 0 as 


Ge (0) — Qn , Gjvf (0) — Qn + i^N = Pn , 


where Pn is the total magnetic moment of the nucleon. Isoscalar (isospin / = 0) and isovector (isospin / = 1) 
components are obtained by the following combinations of proton and neutron FFs 


.s Fl2 + Fi,2 

^ 1 , 2 = -p;- 


_ ^E,M + Ge,M 
^E,M — -;;- 


jpP jpn 

.V ^ 1,2 --^ 1,2 


FE = ■ 


_ '^E,M '^E,M 

^E,M — - 5 - 


(5) 


1.2. Time-like form factors 


In case of annihilation, fig. 1 horizontal direction, following the notation of eq. (1), the amplitude for the nucleon- 
antinucleon production, 7 * —^ NN, is 


(IV(p)iV(p')b'^(0)|0) 


f^N, 2n /j , .fr qv , 2n \ / 

-F, (g) u(p), 


«(p) Fi {q + * 


2Mn 
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where Fl^(q^) and F 2 {q^) are the Dirac and Pauli FFs in the TL region, as indicated by the over-bar. Even in 
this case, the hermiticity of and the time reversal symmetry would imply real TLFFs. However this would be 
true only if the TL photon had not enough virtual mass, q^, to produce physical particles as intermediate states. 
Otherwise, when the values of exceed the production threshold of the lightest allowed intermediate state, the 
amplitude, and hence the TLFFs, become complex. The rising of a finite imaginary part is a consequence of 
unitarity and can be formally demonstrated by considering the optical theorem. 

Since the lightest hadronic physical (on-shell particles) state, allowed by quantum number conservation, is the two- 
pion one, the imaginary part of the amplitude is different from zero starting from the so-called theoretical threshold 
9theo = where M,r is the pion mass. In light of this non-vanishing imaginary part, the hermiticity of the 

current operator and the time reversal symmetry enforce, for the FFs, instead of reality, the Schwarz reflection 
principle and hence a discontinuity across the half line (?theo! 00 ). Such a portion of the TL region is then excluded 
from the analyticity domain or, in other words, it represents a branch cut. 

From the experimental point of view, the extraction of TLFF data involves additional difficulties with respect to 
SLFFs. First of all, TLFFs are complex so, to have a complete determination, moduli and phases, or imaginary 
and real parts, should be measured. However, even by using polarization observables [25], only relative phases, 
between electric and magnetic Sachs FFs, are accessible. 

Moreover, since TL data are extracted from the cross section of the annihilation processes e''"e“ NN, TLFFs 
can be measured only for values above the so-called physical threshold q^^ys = {2 Mn)^- It follow that the TL 
interval [0, (2M]v)^], where TLFFs are still well defined and receive also the most important contributions from 
hadronic intermediate states, is not experimentally accessible and for that reason it is called ’’unphysical region”. 
As already stated [26], taking advantage from crossing relations, SLFFs and TLFFs are interpreted as limit values, 
over the negative and positive real axis, respectively, of unique functions, Fi^ 2 {q^), defined in the whole g^-complex 
plane with the discontinuity cut (^theoiOo), due to unitarity (optical theorem). In more detail 

Fi^iq^) = lim F^(q^ ± ie) q^ < 0 (SL region) 

e->0 

< Fi^{q^) = lim fI^{ q^ ±ie) Q < q^ < ql^^^ (TL region) , * = 1,2, 

Fl^iq^) = lim F^{q^ + ie) q^ > gt\eo (TL region) 

E—>0 


where: the limits in the SL region and in the portion of TL region up to the theoretical threshold can be taken 
indifferently from above or below the real axis, because there is no discontinuity there, while the limit values of FFs 
around the cut, in the TL region, depend on which edges of the cut is considered. In particular, as a consequence 
of the Schwarz reflection principle, the FF values in the upper and lower edges are complex conjugates, i.e., as 
e ^ 0+, 


Fi, 2 {q‘^ + ie) = F’i* 2 (<J^ - ie ), q^ > • 

1.3. Analytic properties of form factors 

Analyticity and unitarity, as well as perturbative QCD, determine important model-independent features of FFs 
(some of which have already been touched upon in the previous section). Any reliable model of FFs must be able 
to reproduce such fundamental features, that concern the analytic structure of FFs as functions of the complex 
four-momentum square and their asymptotic behavior, i.e., the power law that rules their vanishing as |g^| —^ 00 . 
Below we list, without proof, the main properties of FFs that will be addressed and discussed in the next sections. 

• Form factors are function of (f, analytic in the whole complex plane except for the branch cut (qtheo,oo). 
Physical FFs are defined as the values of such functions for real q^. Moreover, from the experimental point 
of view, FFs are measurable only for < 0 (SL region), and q^ > gphys (a subset of TL region). The TL 
interval (0, gpj,yg), being experimentally forbidden for FFs, is called unphysical region. 
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• At high momentum transfer we can invoke the perturbative QCD or the quark counting rule [28] to infer the 
FF asymptotic behavior. In particular, in the scattering channel in order to maintain the nucleon entirety, 
the four-momentum transferred by the virtual photon must be shared among the three valence quarks via 
gluon-exchanges. The minimal number of gluons to be exchanged is two and hence the FFs must contain 
terms with, at least, two gluon propagators that entail the power law behavior 


^E,M 


{q^)~ 



2 , 

q —>■ —oo , 


where the limit is in the SL region. However, such a power law can be extended also to the TL region 
by considering the Phragmen-Lindelof theorem [29], that applies to FFs because of their analyticity and 
boundedness. 

• A very powerful consequence of analytic properties of FFs is the possibility of using a particular analytic 
continuation tool based on the Cauchy theorem [30], i.e., the dispersion relations for the imaginary part 



Im [F{q'^)\ 

ql2 _ g2 _ jg 


dq'‘ 


^ (g?heo, oo) , 


where the symbol F stands for a generic FF. The threshold value to be used as lower limit of the dispersion 
relation integral depends on the isospin of the considered FF. In case of isovector components, intermediate 
states with only even numbers of pions are allowed, hence, as already seen, the threshold is ?theo = (2M^)^, 
while for isoscalar components q[heo = However, using the lower threshold is always correct, since 

the imaginary parts of the isoscalar FFs are vanishing for q^ < q'theo- 

• Another interesting issue, that emerges by considering the definition of and G^,j, and assuming analyt¬ 
icity for the Dirac and Pauli FFs, is the identity G^{AM%) = Gm{4M^). On the other hand, the electric 
and magnetic FFs could be different at the physical threshold only if Fi' and Fi' were singular there [20]. 
Such an identity implies that, at the threshold Qphysi the nucleon vertex is described by a unique FF, i.e., 
there is only one degree of freedom and the cross section, loosing its dependence on the scattering angle, 
becomes isotropic. In other words, even though angular momentum conservation allows S and D waves for 
the NN system produced by one virtual photon (Born approximation), at the production threshold the 
D-wave contribution must vanish, so that only the isotropic S-wave survives. 

In principle, the identity G^{4M%) = Gm{4M^) can be verified experimentally by measuring, for instance, 
the ratio Ge /Gm at the physical threshold. However, in a symmetric e''"e“ collider, the possibility to reach 
or, even, get very close to the threshold, is prevented by physical limitations. Indeed, in case of the annihi¬ 
lation process e''"e“ —>■ pp, the proton and the antiproton, produced almost at rest in the laboratory frame 
would have no enough momeutum to reach the detector. 

In the last fifteen years, the so-called ’’initial state radiation technique”, developed at the flavor factories, 
allowed to avoid this limitation, so that values of the ratio Ge/Gm [31] have been measured very close to 
the physical threshold. These data, together with older measurements performed in the crossed channel 
pp —>■ e'''e“ [32], agree with threshold-isotropy requirement G^{4M^) = G\^{4M%), but do not exclude 
possible, small D-wave contributions. 


2. Skyrme model 

The Skyrme model was introduced by Tony Skyrme in 1960 as a model for strong interactions [19]. The basic and 
innovative idea was that fermions could emerge as particular, stationary and quantized solutions of a non-linear 
field theory with only boson fields. Stationary solutions of this kind are usually called solitons, the quantized 
ones, associated to the Skyrme Lagrangian, are instead called skyrmions. 

The interest in this model increased when’t Hooft and Witten proposed the l/N^ expansion of QCD [33, 34] and 
Witten showed that the Skyrme model led to a Lagrangian which was equivalent to that of the 1/A^c expansion. 
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The first application of this model is due to Adkins, Nappi and Witten [35, 36], who computed some static 
quantities for nucleons by obtaining a quite acceptable (~ 30%) agreement with the measured values. Such an 
agreement strengthened the conviction to being on the right track to achieve an effective approximation of QCD 
at low energy. 

In order to build up a representation of nucleon SLFFs, we will follow the work of Braaten, Tse and Willcox [21], 
that, in 1986, for the first time, used the Skyrme model to compute nucleon FFs. 

2.1. Skyrme Lagrangian 

The Skyrme Lagrangian, which is based on the Lagrangian of the so-called non-linear a-model [37], has an 
SU{2)l X SU{2)ii chiral symmetry which is spontaneously broken to SU{2). Assuming that the isoscalar and 
isovector fields a and if, as a consequence of the symmetry breaking, are linked by the relation 


where FV = 108 MeV ]36] is the weak pion decay constant, see tab. 1, the Skyrme Lagrangian can be written in 
terms of the only SU{2) field 


U{r) = ^ ^o'(r) -\-iT ■ 7f(a:)^ = exp (ir ■ F{f)^ 


( 6 ) 


where: f is the vector of Pauli matrices and the function F{r) is the ’’axis-angle” representation of the chiral field 
U{r). The complete Skyrme Lagrangian density, that will be used in the following, reads 


-Cskyrme = ^ Tr17+) + ^ [[D^UU\D,UUY) + Cwz + Tr (l7 -f t/+ - 2j 


(7) 


is the covariant derivative, which includes the electromagnetic interaction. Besides the usual kinetic term, 
the second contribution, which is quadratic in the field derivative and represents a repulsive short-range potential 
with a coupling g, has been introduced ad hoc by Skyrme in order to have stationary solutions. The contribution 
Cwz, called Wess-Zumino term [38], which accounts for the QCD anomalies, is written as a non-gauge invariant 
coupling between the photon and a conserved topological current [39], i.e., 

Cwz = , with: S'" = ^e^"^'"Tr(f/+a.f/aAl7+a,,17). (8) 


The topological charge associated to corresponds to the baryon number B, hence the baryons are identified 
as those solutions with B = 1. Finally, the last contribution of eq. (7) is a mass term, which explicitly breaks the 
chiral symmetry and it is treated perturbatively. 

We consider the particular class of solutions obtained by specializing the axis-angle function, F{r) of eq. (6), 
according to the hedgehog ansatz: F{r) = F{r)r/\r\ = F{r)r. In this way, the space of the surviving 517(2) 
symmetry, which is the isospin, takes the radial configuration of the r-space. In other words, in a given position 
r, the isospin vector has the same direction and orientation of the position vector r. The intensity of the axis- 
angle function F(r), indicated with the symbol F{r), is called chiral angle. Stable solutions, the skyrmions, 
are stationary minima of the energy, obtained by solving the Euler-Lagrange equation, which is a non-linear 
differential equation for the chiral angle F(r). 

The further step, that is the skyrmion quantization, consists in quantizing collective modes, translations and 
rotations, in the isospin space. This can be done by using, in the Lagrangian density Bskyrme, instead of the field 
U{r) of eq. (6), its time-dependent version 


U{r,t) = Ait)U{r-X{t))AHt), 
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where A{t) is a uniform SU{2) matrix and X{t) is the skyrmion center-of-mass position vector. As a consequence 
of rotational and translational invariance, the resulting Lagrangian depends only on the derivatives of A{t) and 
X{t), and it reads 


iskyrme — J T -Cskyrme — —M H-^-h ATr(j4^A) , (9) 

where M and A are mass and moment of inertia of the skyrmion. 

Owing the hedgehog ansatz, the rotational operator is related, not only to isospin, but also to spin. The skyrmion 
can be interpreted as a nucleon by requiring the rotational operator to have a semi-integer eigenvalue, so that, 
spin and isospin are both quantized to 1/2 [40]. The Hamiltonian for the quantized skyrmion, in terms of its 
three-momentum and spin operators P and S, is 

Hs^.yrme = M+^P^+^^. ( 10 ) 

Since, P, S 3 and I 3 (third components of the spin and isospin) are mutually commuting operators (they also 
commute with the Hamiltonian), the system described by the Hamiltonian of eq. (10) is manifestly non-relativistic 
and has the eigenstate |p, 53 , 13 ), where p, S 3 and 13 are the corresponding eigenvalues. Moreover, as already 
noticed, there is also the conserved topological charge B. Thus the nucleon is identified as the eigenstate of 
Hskyrme with: S 3 = ts = 1/2 and 5 = 1. 

The skyrmion mass M, appearing in eqs. (9) and (10), represents the minimum of the energy, obtained by solving 
the Euler-Lagrange equation of the Skyrme Lagrangian, eq. (7), to find the solitonic solution, i.e., the chiral angle 
Fir). 


Quantity 

Unit 

This work 

Experimental 


MeV 

(fixed) 

108 

~ 186 


MeV 

(fixed) 

138 

~ 138 

g 


(fixed) 

4.84 

- 

M 

MeV 

937 

~ 938 


fm 

0.88 

0.8775(51) 


fm^ 

-0.31 

-0.1161(22) 

1/2 

V P/M 

fm 

0.79 

0.777(16) 

\' n/M 

fm 

0.82 

0.862(9) 

Pp 

/is 

1.97 

2.792847356(23) 

/In 

/is 

-1.24 

-1.9130427(5) 

/Ip / f-^n 


-1.59 

-1.459898075(5) 


Table 1. Parameters and static quantities obtained in the framework of the Skyrme model compared with their experimental 
values [41]. Our results, being obtained in the same conditions, reproduce quite well those of Ref. [36]. The ~l%o 
difference in the nucleon mass is probably due to a slightly different normalization range for the chiral angle F{r). We 
used 0 < r < 8 fm, while, in Ref. [36], there is no indication on that. 


2.2. Electromagnetic form factors in Skyrme model 

Nucleons FFs were firstly obtained in the framework of Skyrme model, by Braaten, Tse and Willcox [21]. The 
starting point consisted in deducing the most general expression for the electromagnetic current, at a given order 
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in some expansion parameter, that fulfilled all symmetries and constraints of the model. Due to the equivalence, 
discussed in Sec. 2, of the Skyrme and an SU{Nc)-QCT) Lagrangian, in the large-limit, the natural expansion 
parameter turns out to be l/N^. For instance, the Skyrme Hamiltonian of eq. (10), being both M and A of order 
Nc, contains terms up to the first order in the 1 /Nc expansion. 

The expression of the electromagnetic current, at the leading 1/Nc order, is written in terms of position, mo¬ 
mentum, spin and isospin operators, and also two model-dependent four-vector functions of (eight scalar 
functions), to be furthermore specified by imposing symmetries (hedgehog ansatz included) and physical con¬ 
straints (e.g.: B = 1). Only two, b{X^) and t{X^), out of the eight scalar functions, survive the characterization 
procedure, hence, using relations and definitions of eqs. (3) and (4) in the Breit frame, nucleon SLFFs can be 
obtained as the Fourier transforms: 


G|(Q^) = y d^rjo{Qr)b{r), 

(11a) 

Gb(Q^) = ^ j if rjo{Qr)r^t{r), 

(11b) 

Gm(Q^) = ^ J f rji{Qr)rb{r), 

(11c) 

Gm(Q^) = |^ j f rji{Qr)rt{r), 

(lid) 


where ji{x) is the i-th spherical Bessel function and = —(f, with = Q > 0 in the SL region. 

In the Skyrme model, the functions b(r) and t{r), are defined only in terms of chiral angle F{r). Their expressions, 
obtained by comparing the electromagnetic current, as extracted from the Lagrangian of eq. (7), with the general 
"educated” parametrization discussed so far, are 


F'{r) sin^ [F’(r)] 

2tx^ ’ 

^ F^ sin^ [F(r)] ^ 1 sin^ [Fjr)] / xi 2 sin^ [F{r)] \ 

4 J .2 g 2 ^2 ^2 j- 


( 12 ) 


The chiral angle F{r) is obtained as solution of the Euler-Lagrange equation, a non-linear differential equation 
of second order, which follows from the functional minimization of the skyrmion mass-energy computed with the 
Lagrangian density Tskyrme of eq. (7) and the representation of the chiral field U{r) given in eq. (6). 

Since the function F{r) and hence b{r) and t{r), are known only numerically, the Fourier transforms of eq. (11) 
allow to compute nucleon FFs only in their convergence domain, which corresponds to the SL region, i.e., > 0 

and Q > 0. Indeed, negative values of would imply pure imaginary \/C^, so that the Bessel spherical functions 
become exponentially divergent as r —>■ oo. 

It follows that there is no possibility to perform analytic continuations of the SLFFs in the TL region, if b{r) and 
t{r) are known only numerically. In Sec. 3 we will develop a procedure to overcome this limitation. 

The order in the l/N^ expansion of the FF expressions given in eq. (11) can be easily inferred by the presence of 
factors M and A at numerator or denominator. In particular, the isoscalar electric and magnetic FFs, eqs. (11a) 
and (11c) are of order zero, the isovector electric FF, eq. (11b) is of order one, while the isovector magnetic FF, 
eq. (lid), is of order minus one, i.e. 0{Nc). Such a heterogeneity should prevent the possibility to combine these 
expressions, following eq. (5), to obtain proton and neutron FFs, as actually has been done in Sec. 3.3, in order to 
compare our results with the available data. The fair agreement, which has been obtained, could be an indication 
that the non-leading contributions to the nucleon FFs, in the 1/A^c expansion, are sub-dominant as if there were 
additional suppression factors. 


2.3. Relativistic corrections 

To extend the results obtained for the SLFFs to high , relativistic corrections have to be included. However, 
the procedure to obtain relativistic skyrmions is still debated and different methods are present in literature. To 
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include relativistic corrections in our FF parametrizations we will follow Ref. [42], In particular, in the SL region, 
relativistic FF expressions are obtained from the non-relativistic ones as 


\ 1 + 


^]V,rel/^2N _ 1 / Q 

VV J — - Lta 


1 + 


1 + 


(13) 


while, in the TL region [20], using = —Q^, 




/ 2 
GEiq 


GT\q^)={ 


iMi 


GUq") 


-2 


q^ < AM% 


q^ > AMh 

q^ < AM% 


q2 


amI 


AMl 


--2 


q^ > AMjf 


(14) 


where the relativistic forms are labelled by the superscript ’’rel”. 

It is important to stress the fact that, unless one considers fine tuned, non-relativistic FF expressions, with zeros 
of particular orders at particular finite values of Q^, such corrections are incompatible with perturbative QCD 
predictions concerning the asymptotic behavior. 

Moreover, since the asymptotic SL and TL limits for a given FF are different, for instance G^’^®*((5^) tends to 
the non relativistic value G%{AM%) as ^ oo and to G^{%M%) as —>■ —oo, these corrections do not even 
verify the Phragmen-Lindelof theorem. 

The effect of relativistic corrections is shown in £g. 2, in case of SL electric and magnetic proton FFs. Their 
inclusion improves the agreement with data at higher , even though, at very high momenta, as expected, the 
agreement worsens. 


3. Analytic continuation and results 

Once the profiles b(r) and t(r) are known, FFs are computed as their Fourier transforms. However, since only 
numerical expressions of b(r) and t(r) can be obtained, the applicability of the representations given in eqs. (Ha¬ 
lid) is limited to their domain of convergence, i.e., the SL region. Indeed, TL transferred momenta would imply 
divergent real exponentials in the Fourier integrals. The only way to overcome such a limitation and hence to ob¬ 
tain FF values also in the TL region, consists in performing an analytic continuation of the original representation. 
The simplest procedure is represented by a direct, analytical computation of the integral, which would return, 
for the FFs, well defined expressions, depending on the variable Q^, that are manifestly analytic. In this case the 
direct approach is prevented by the lack of analytic forms for the profiles b(r) and t(r). Nevertheless, the fact that 
they are well known in a wide range of r, allows us to define simple analytic functions that approximate them 
with an accuracy that, in principle, could be indefinitely improved by increasing the number of free parameters 
to be settled. Moreover, the structure of such fit functions is inferred by the knowledge of the profiles in the two 
limits: r —> 0 and r —> oo. To expound the analytic continuation procedure that we have developed, we consider 
in detail the case of the electric isoscalar FF G%{Q^). Its integral representation, given in eq. (Ha), can be also 
written as 


G%{Q^) = ^ ^ ® rb{r)dr 


(15) 


and it contains only the profile b{r). 

The behaviors in the limits r —>■ 0 and r —>■ oo are well known because the corresponding differential equations 
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(GeV^) (GeV^) 

Figlir© 2. Left panel: proton electric FF. Right panel: proton magnetic FF normalized to the proton magnetic moment. Solid 
and dashed lines represent the numerical results with and without relativistic corrections, respectively, while the empty 
circles are the world data sets [7]. 


can be analytically solved and we get 

b(r) ~ /lo + h2r‘^ , 

r -»-0 


b{r) ~ /icx) 

r—¥oo 


-3M^r 1 


(16) 


where ho, ho and hoa are free parameters. The exponential, which guarantees the fast vanishing of the profile, 
plays a crucial role in characterizing the analytic structnre of the FF, hence the fit function has been defined as 


foflt(r) 


Pn{r) -3M„r _ -3M„r 

Pr.{r)^ ^Y.%ohri^ 


(17) 


where Pn{r) and Pm{r) are polynomials of order n and m respectively (m, n G N), with real coefficients 
and {bj}JLo, so that: a„ yf 0, 6m 7^ 0 and bo yf 0. 

The conditions given in eq. (16) imply: (oi — aobi/bo — 'iM-nao/bo) = 0 and m — n = 5, respectively. The rational 
part of bfit(r) is a meromorphic function with a finite number M < m oi distinct poles in the r complex plane, 
hence it can be written as the Mittag-Leffler expansion [43] 


Pr.(r) 

Pm{r) 


E 

1=1 




-flj 


where {zj}fLi C C and C N represent the poles and the corresponding multiplicities, while is 

the fc-th coefficient of the Laurent series about the pole Zj. Assuming that Pm{z) has only order-one zeros not 
coincident with those of Pn{z), the function bfit(r) has only simple poles, i.e., /ij = 1 for all j = 1, 2,, M, and 
M = m. In this case the Mittag-Leffler expansion reduces to 




(18) 
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where the coefficient Rz^ is the residue of the j-th pole 


Rz^ = Ci'] = Res 




P™(r) 




E n 

i=0 


T.i:=,kb^^T 


j = l,2,...,m. 


(19) 


We have an additional condition on the parameters zj: they can not be positive real numbers, because b(r) has 
no poles for r > 0. 

The fit function 6flt(r), per se, has apparently no physical content, because its parameters are not directly 
connected to physical properties of the system under consideration. Nevertheless, as will be discussed in the 
following, the poles Zj play a crucial role, indeed they define the analytic structure of FFs in the complex 
plane. 

The orders m and n of the polynomials and hence the number of free parameters that define 6flt(r), have been 
chosen by following a criterion which combines the higher accuracy in the description with the smaller redundancy 
of parameters (a given parameter is defined redundant when its inclusion does not improve the accuracy of the fit). 
The only constraint on the orders m and n of the polynomials is about their difference, that must be: m — n = 5. 
A satisfactory ht has been obtained with n = 3 and m = 8. The resulting electric isoscalar FF is 


gUq") = ^ 


2'KRzj 

iQ 


■ roo (iQ — 3MTT)r 

roo {—iQ—3MTr)r 

/ ® dr 

/ ® dr 

1 

o 

1 

O 


( 20 ) 


where, having only poles of order one, 

Rzi = Res 


rPn{r) 


Pm{r) ’ 

All integrals appearing in eq. (20) belong to the same class 


H{ap) = 


f 


r + /3 


dr , with: 


— Pzj • 


Re(a) > 0 

P ^ (-00,0] 


( 21 ) 


( 22 ) 


The conditions on the parameters a and /? ensure the convergence of the integral that, as can be easily seen by 
making the substitution w = ar, depends only on the product a/3. In particular, the integrals of eq. (20) can be 
obtained with: a = 3M,r ± iQ and /3 = —Zj. These assignments, having no poles on the positive real axis and 
being Q > iM-jr, automatically fulfill the convergence conditions. 

In the a/3 domain defined in eq. (22), the function H{aP) has also the following representation (see App. A) 


H{aP) = e°‘^Ei{aP) = e°‘^ 


-7 - ln(a/3) + ^ 

k=l 


(-!)'=+! (a/S)* 
kk\ 


(23) 


where E-i_{z) is the exponential integral function or ” Explntegral” and 7 is the Euler-Mascheroni constant [44]. 
Finally, using eqs. (22) and (23), the representation of eq. (15) can be integrated to obtain 


G|(g) = 5 [(.Q _ [{-iQ - 3M^)zj]'j . (24) 

Since the function Ei{z) is analytic in the whole 2: complex plane with a cut along the negative real axis^, it is 
now possible to extend the parametrization for G% to the TL region, by making the substitution Q —^ iq (q > 0), 
so that 

gUQ = iq) = — -Y. Y K-q - 3M^)zj] - ](q - 3M.)2^]} . (25) 


^ This is a typical logarithmic branch cut as can be seen in the representation of Ei{z) given in eq. (23). 
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3.1. The branch cut in the complex plane 

The properties of the representation obtained for G% and, in particular, the presence of branch cuts, as well as 
their location in the complex plane, depend on the analytic structure of the Expintegral functions. Following 
the derivation given in App. B, we obtain for G% in the SL region the expression 

G%{Q) = Rr,lm{H [{iQ - 3M^)rj]} + ^ ^ Im [(iQ - SM^)cj] - Rc,H[{iQ - 3M^)*cA (26) 

^ .--1 ^ j^l ^ ^ 

e(Q- 3M, d{yj)R,^ g(iQ-3M,)c, ^ e{-yj)Ro, , 

where rj and Cj are real and complex poles, Xj and yj, in the arguments of the Heaviside theta functions, represent 
real and imaginary parts of Cj, while Rr^ and Rc- are the residues. For a detailed description see App. B. The 
parametrization of eq. (26) is explicitly real in the SL region, i.e., for Q > 0. In the TL region the G% expression 
becomes 




GE{iq) 


—^ ^ Re {Rc^H[{-q - 3M^)cj] - R,^H[{q - 3M,)c,l} 
t=i 
2 ' 

-I] A, {H [(-q - 3M^)rj] -H[{q- 3M^)rj]} 

t=i 

+ — - 3M^)e{-Xi) [Rc,e^'^-^^-'>^^9{yj) + 



(27) 


While the first term is real, the second and the third could have non vanishing imaginary parts. In particular, 
the second term contains H functions, that embed a logarithmic structure and hence, for negative arguments, 
have non-zero imaginary parts. Having only negative real poles, Vj < 0, and q> 0, the argument (—g —3M^)rj is 
always positive, whereas (g —3M,r)rj can be negative when q > 3M.^, following exactly the theoretical requirement 
(see Sec. 1.3). A further imaginary contribution, given by the last term, is due to non real poles with negative 
real part. Finally, the TL imaginary part of G%, which is non vanishing only if there are poles with negative real 
parts, is given by 


Im[G|(ig)] 


- e{q - 3M.) 

<7 


^6i(-a;j)Re 


0=1 


LRc 


(g-3M7r)c 


1 ^ 

loE 


.7 = 1 


Rrj e 


(q-3M^) 



In summary: 

• the profile b{r) is obtained as numerical solution of a differential equation; 

• such a solution is fitted with bnt(r), a product of a rational function and an exponential, that fulfills the 
requirements for r —>■ 0 and r —>■ oo and moreover it has only simple poles not belonging to the positive real 
axis; 

• the rational part of bfit(r), being a meromorphic function, can be written as a Mittag-Leffier sum and hence 
its Fourier transform, which is the FF G%{Q), is a sum of Fourier transforms of simple poles, Zj, multiplied 
by an exponential, i.e., Expintegral functions with arguments: (±iQ — 3M,r)«j; 

• the analytic properties of G%{Q), especially when there is at least one pole with negative real part, reproduce 
exactly those expected for FFs on the basis of first principles. 
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The isovector FFs are obtained with the same procedure described in detail for G%, i.e., by fitting the profile func¬ 
tion t{r) with a ratio of polynomials and an exponential deduced from the solution of the asymptotic differential 
equations. However, in this case, to account for the two and four-pion coupling, two exponential contributions 
are considered 


t{r) 


P„(r) -2M^r . Pn{f) -4M„r 

PArf PUrf 


Following the line of reasoning used to study G%, this expression leads to two different branch cuts, that generate 
from the two thresholds: Qtheo = (2M^)^ and g^theo = 


3.2. Asymptotic behavior 

As already discussed, the asymptotic behavior of the FFs obtained with this procedure is completely determined 
by the relativistic corrections of eqs. (13) and (14). Nevertheless, it is interesting to study the high-Q behavior of 
the non-relativistic (uncorrected) FF expressions. Such a behavior can be derived from the asymptotic expansion 
of the Expintegral function [45] 


= V J +0[{n- l)!|zr"] , 2 ^ oo. (28) 

fc =0 

The rigorous treatment is given in App. C, where the SL and TL asymptotic behaviors are obtained for general 
profile functions, but not taking into account the branch cut corrections discussed in Sec. 3.1. However, as we 
will see in detail in the case of G%, such corrections do not spoil the power law behavior driven by the expansion 
of eq. (28). Using eq. (C.13), the SL isoscalar FF in the high-Q regime can be written as the series of increasing 
powers of Q”'^ 


OO 

G%{Q)- EssnW)- 

Q—^oo ‘ ^ 
fc=0 


(29) 


with: (?g'2(Q) ~ Q ^ and g^gliQ) ~ (Q ^^2int[{fc-i-i)/2] ^ > p jjj particular, when Q —^ oo, the first four terms 

behave as 





247r Mtt 

Q" 


^Re 

i=i 


2887r Mtt 


^Re 

1=1 



(30) 


Each of them depends on the corresponding derivative (the fc-th derivative for the fc-th term) of the rational 
function r evaluated in the origin, i.e.. 


AL 

dP 




d’^ rP„(r) _ d’‘ ^ Rj 
dr'^ Pm (r) ^ i* - Zj 



1=1 



fc = 0,1,... . 


Since the function r feflt(r)e®*^’'’' vanishes in the origin, having no poles there, the first contribution {k = 0) in 
eq. (30) is also vanishing and hence g^sliQ) g^sliQ)^ ^^e same order in Q, i.e. Q~‘^, become 

the leading terms. The TL expression for G%{iq), given in eq. (27), apart from the factor q~^, has two kinds 
of contributions: the first depends on the functions H{z), while the second, which accounts for the branch 
cut corrections, has an exponential behavior. More in detail, there are two exponentials that, being complex 
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conjugates, have the same modulus and hence the same asymptotic behavior. Their moduli scale like ~ 
when q —>■ oo, where xj is the real part of the j-th pole. However, such contributions are weighted by three 
Heaviside theta functions, one of which ensures the strict positivity of Xj , hence all the exponentials are vanishing 
as q —^ oo and the asymptotic behavior is dominated by the only terms which contain the H{z) functions. In 
light of that, the asymptotic behavior of G%{iq) can be described in terms of the series 


OO 

G%{iq) ~ 

o—>-oo ' ^ 

fc =0 

similar to that of eq. (29), where the functions g^liq) are defined by the TL expansion of eq. (C.15), and are 
related to the corresponding SL terms as g^l{q) = other words, the TL asymptotic behavior follows 

the same power law as in the SL region. The leading contributions, in both regions, are determined by the 
behavior of the profile function in the origin r = 0. 

However, while for the electric isoscalar FF the obtained behavior agrees with the perturbative QCD prediction, 
i.e., G%{Q) ~ Q~‘^ and G%{iq) ~ , for all the other three FFs, see eqs. (C.17) and (C.19), we achieved 

Q—>-oo g—^oo 

the faster vanishing behaviors 

GUQ),Gl,{Q),Glt{Q) ~ Q-\ Gl{iq), Glt{iq), G\t{iq) ~ g"®. (31) 

Q— q—^CG 


The only possibility to recover the expected power laws should be that to consider a profile function having in 
the origin a zero of a lower order. For instance, in case of GXi, the profile function is /(r) = r^t{r), see eq. (lid), 
and, as r —>■ 0, /(r) oc r*, with I = 3, because the function t{r) is finite and non vanishing in the origin. This 
power, I = 3, determines (see eq. (C.18)) the asymptotic behavior as given in eq. (31). On the other hand, the 
perturbative QCD expectation, i.e., the power laws Q~‘^ and q~'^, in SL region and TL region respectively, would 
be obtained only with I = 2, which means that t{r) should have a simple pole in the origin. 

3.3. Results 

To have a direct comparison with data, results are given for the electric and magnetic Sachs FFs of proton and 
neutron, and G’e,m! even though the primary outcomes of this procedure, see eqs. (11), are their isospin 

components G% and G^ These two sets of FFs are related by the linear combinations given in eq. (5). 
Moreover, SLFFs and TLFFs will be given as functions of and q^ respectively, with the simple convention 
G^:#^(±q^) = G|'"#'^(tQ^). 

For each FF both relativistic corrected and uncorrected results are presented as solid and dashed curves, respec¬ 
tively. In the TL region below the physical threshold ^phya, the two curves coincide, see eq. (14). 

Moreover, all the SLFFs, but for G^ which is vanishing at = 0, are normalized to the so-called dipole FF 

Gr,(0^)= (^1+^) , (32) 

with M|) = 0.71 GeV^. Such a FF, with only one free parameter, the dipole mass Md, describes quite well SL 
data on G^{Q^), G^j{Q^)/g-p and G'hiQ'^)/gn, as can be seen in fig. 3 and fig. 5 right panel, where indeed the 
data (empty circles) spread out around the unity. 

3.3.1. Space-like region 

Figures 3-5 show predictions (solid and dashed curves) and data (empty circles) for the SL electric and magnetic 
FFs of proton and neutron. In particular, solid and dashed curves represent the predictions including and not 
including the relativistic correction described in eq. (13). In the case of the proton, fig. 3, the relativistic- 
uncorrected predictions, also thanks to the constrained unitary normalization at — 0, describe quite well data 
up to ~ 0.4 GeV^. Above this limit the predictions start to decrease faster than the dipole. Such a behavior is 
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(GeV^) 



(GeV^) 


Figure 3. The Electric (left panel) and magnetic (right panel) proton FFs in the SL region, normalized to the dipole FF and 

also to the magnetic moment, are compared with the world data sets, empty circles, from Ref. [7]. The solid and dashed 
curves represent the predictions for the FFs obtained including and not including relativistic corrections as given in 
eq. (13). The data sets are the same of fig. 2, however, due to the different scales, logarithmic and linear, and to the 
dipole normalization, the errors appears larger. 




Figure 4. Electric (left panel) and ratio electric to magnetic (right panel) proton FF in the SL region. The ratio is compared with 
the data, empty circles, obtained by means of polarization observables [7]. The solid and dashed curves represent the 
predictions that include and do not include relativistic corrections, respectively. The vertical line (red in the on-line 
version) in the left panel indicates the value — 4Mfj. 


expected in case of the magnetic FF, in fact, as shown in eq. (31), the power law that rules its high-Q^ vanishing 
is On the other hand, the electric FF, due to the contribution of G%, see eq. (C.17), should tend to zero as 
Q“^, i.e., at the same rate as the dipole. 

The obtained faster vanishing behavior is due to the presence of a zero for at Qq ~ 2.3 GeV^, see 

the dashed curve in left panel of fig. 4, so that G^{Q^) —>■ 0“ (from below), as O when > oo. The 

agreement with data is improved by including the relativistic corrections, solid curve in fig. 3. In particular in case 
of G^, left panel of hg. 3, the prediction follows the trend of the data, i.e., the dipole behavior, up to — 2 GeV^ 
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and then it drops down. This is a conseqnence of the Q^-dilation natnre of the relativistic correction, that moves 
the zero for G^{Q^) from Qq to Q^ei — 4M^/(4M^/Qo — 1) — 8.4 GeV^, see the solid curve in the right panel 
of hg. 4, and hence the quick descent is shifted at higher Q^. 

As already discussed in Sec. 2.3, the asymptotic behavior of the electric FF is drastically modihed by the relativistic 
corrections, in fact, tends to the finite value G^(4M^), i.e., 

^ G|(4M^) ~-0.0025. 


The fact that such a value is very close to zero, see the vertical line in the left panel of hg. 4, and that the 
uncorrected electric FF scales as the dipole makes the corrected FF closer to the data. 

Also the prediction for the magnetic proton FF, right panel of hg. 3, improves its agreement with data up to 
~ 3 GeV^, when the relativistic corrections are considered. In this case, at high , the prediction gets larger 
than data, see also the right panel of hg. 2, and its steep rising, from ~ 4 GeV^, is a consequence of the 
normalization to the dipole. Moreover, asymptotically goes like Q~^, so that the ratio to the dipole grows 

like Q^. Contrary to the case of G^, no zeros are found for G^^. 

The right panel of hg. 4 shows the ratio between electric and magnetic proton FFs normalized to the proton 
magnetic moment, the solid and dashed curves are the predictions with and without relativistic corrections, 
while the empty circles represent the data extracted from polarization transfer observables in e-p scattering [7]. 
Such experimental values show an unexpected linear decreasing trend, whose extrapolation would give a zero at 
~ 10 GeV^, which is close to the obtained value Qo.rei — 8-4 GeV^. 




Figure 5. Electric (left panel) and magnetic (right panel) neutron FF in the SL region compared with the worid data set, empty 
circles, from Ref. [7]. The solid an dashed curves represent the predictions for the FFs obtained including and not 
including relativistic corrections as given in eq. (13). 


Electric and magnetic FFs of neutron are shown in fig. 5 in comparison with the data. The two predictions, also 
in this case, refer to the relativistic corrected (solid curve) and uncorrected (dashed curve) results. Apart from 
the low-Q^ region, where the normalization forces the predictions to follow the experimental points, the agree¬ 
ment with data appears worse with respect to what has been found for the proton. The inclusion of relativistic 
corrections does not improve the accordance with data, in particular, the agreement is even worsened in case of 
GJ, left panel of hg. 5. Finally, no zeros are found for G^ and GJ/. 

Concerning the asymptotic behavior of neutron SLFFs, the same conclusions driven for the proton can be consid¬ 
ered. In particular, as —>■ oo, the uncorrected predictions for G^ and G^i scale as Q~'^ and Q~^, respectively. 
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while the corrected behaviors are 


r GY^\Q^) 
\ G:r\Q") 


G'S(4M|,) ~ 0.012 


Q 


2 G'm(4Mjv) 


3.3.2. Time-like region 

Results and data in the TL region will be described as functions of the positive, squared four-momentum transfer 
= —Q^ > 0. As extensively discussed in Sec. 1.2, starting from the theoretical threshold f/theo = (2M^)^, 
FFs develop non-vanishing imaginary parts due to the coupling of the virtual photon, which now has enough 
virtual mass, with hadronic intermediate states. It follows that, in this kinematical region, the nucleon structure 
is described by four real functions, i.e., real and imaginary parts of the electric and magnetic FFs. 



(GeV^) (GeV^) (GeV^) (QeV^) 

Figure 6. Real and imaginary parts of the proton electric FF (a,b) and magnetic FF, normalized to the magnetic moment, (c,d) 
in the TL region. The vertical (red in the on-line version) lines indicate the physical threshold . 





Figure 7. Real and imaginary parts of the neutron electric FF (a,b) and magnetic FF, normalized to the magnetic moment, (c,d), 
in the TL region. The vertical (red in the on-line version) lines indicate the physical threshold — (4Mjy)^. 


Figures 6 and 7 show real and imaginary parts of TLFFs for proton and neutron respectively, including (solid 
curve) and not including (dashed curve) the relativistic corrections, as given in eq. (14). Such corrections become 
effective only above the physical threshold = {2Mn)‘^- For all these quantities there are no available data 
and moreover, even in case of an ideal experiment able to exploit also polarization observables in annihilation 
processes, only relative phases between and Gm would be accessible, besides their moduli. In fig. 8 the 
relativistic corrected, TL (solid curve) and SL (dash-dot curve) moduli of the four nucleon FFs are represented as 
functions of |g^|. Apart from the very first portion of the unphysical region, 0 < < 1 GeV^, where the opening 

of the logarithmic branch cuts manifests itself in bumpy behaviors, TLFFs are smooth decreasing functions of 
q^. Moreover, as it is shown in hg. 8, TLFFs are systematically larger than their SL counterparts at \Q^\ = \q^\. 
Such a discrepancy contrasts with the Phragmen-Lindelof theorem (see Sec. 1.3), stating that SL and TL limits 
of a given FF should correspond. However, on the one hand, as already discussed, relativistic corrections entail 
important modifications of the asymptotic behavior, and on the other hand, it seems plausible to consider as 
center of mass of the SL-TL symmetry not simply — 0 but rather a TL value, say gcMi lying inside the 
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unphysical region. In light of this we should expected Ge,m{Q^) — \Ge,m + 2^0^)I (with the argument 
we mean SLFF at |Q^| = \q^\) and using, for instance, qcyi = 1 GeV^, the SL-TL discrepancy can be reduced. 



Figure 8. Moduli of TL (solid curve) and SL (dash-dot curve) electric and normalized magnetic FFs of the proton (a,b) and 
neutron (c,d). The symbol |g^| stands for positive (TL) and negative (SL) . Only relativistic corrected values have 
been considered. 



Figure 9. Real and imaginary part of electric (dotted curve) and magnetic (solid curve) FFs of the proton (a,b) and neutron (c,d) 
in the TL region, across the physical threshold — 4AL^, indicated by the vertical line (red in the on-line version). 


Electric and magnetic FFs of proton and neutron are obtained using the combinations, given in eq. (5), of the 
isospin components, which are the Fourier transforms, see eq. (11), of the two functions b(r) and t{r), defined 
in terms of the same chiral angle F{r) through the non-linear differential eq. (12). It follows that, FFs are all 
non-trivially interconnected. Moreover, as given in eq. (4), G% and are also linearly related to the Dirac and 
Pauli FFs, in such a way that, assuming no singularity at the physical threshold ^phys = (2Mjv)^ for and 
F 2 , the electric and magnetic FFs of each nucleon must coincide at such a value. As explained in Sec. 1.3, 
the identity G;e(4M^) = Gm( 4M^) implies (it is a sufficient condition for) isotropy at the production threshold, 
i.e., the differential cross section for e''"e“ —f NN in the e''"e“ center of mass frame. 


dcosO 



(33) 


loses its dependence on the scattering angle 0 as f (gphys)^- This also means that, even though parity 
conservation allows S and D-wave for the NN system, at the production threshold only the S-wave can contribute. 
So that, by reversing the argument, the violation of the identity® G%{4:M%) = Gm{4M^) would imply anisotropy, 
i.e., the presence of a D-wave contribution also at threshold or, equivalently, the presence of singularities in the 
Born amplitude. 

Figure 9 shows a comparison between real and imaginary parts of electric and magnetic FFs for proton and 


® Being TLFFs complex functions of the equality G^(4M^) = G((^(4M^) is equivalent to two independent identities 
for the real and the imaginary parts. 
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neutron, in the region of across the physical threshold Qphya (vertical line). To verify the equality = 

the two pairs of real parts (Re[G|(4M|r)] = Re[G^p(4M^)] and Re[GS(4M|r)] = Re[G)(f (4M|r)]), as 
well as the two of pairs of imaginary parts (Im[G^(4M|,)] = Im[G)(f (4M^)] and Im[G£(4M^)] = Im[G)(f (4M^)]), 
would coincide at the threshold. Since no one of these identities is verihed, there is no coincidence between electric 
and magnetic FFs at the production threshold (see next section for a detailed discussion). 




Figure 10. Moduli of the electric (dotted line) and magnetic (solid line) proton (left panel) and neutron (right panel) FFs, in the 
TL region, across the physical threshold — 4M^, vertical line (red in the on-line version). 


It is interesting to notice that the differences among real and imaginary parts at the threshold are partially 
compensated when moduli are taken into account, as shown in fig. 10. Nevertheless, there is isotropy-violation at 
the threshold = gphys a-s it is shown in hg. 11, where, in the left and central panel, are reported moduli of the 
S-wave and D-wave, proton and neutron FFs which are defined in terms of Sachs FFs as 

2y/Wm^)GUe) + GUq^) y/¥J¥^GU<f)-GU<f) 

[Q ) — - o - ’ ) — - n - • 





Figlir© 11. Moduli of the S-wave (dotted line) and D-wave (solid line) proton (left panel) and neutron (central panel) FFs, and 
their ratios (right panel), the solid curve for the proton and dot-dash curve for the neutron. The vertical line (red in 
the on-line version) indicates the physical threshold <?phys ~ 


The right panel of fig. 11 shows the relative contribution, in modulus, of the D-wave with respect to the S-wave 
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FF. It turns out that, in case of the neutron (dash-dot curve), the isotropy-violation is stronger, indeed the D-wave 
dominates in the region around the threshold (/phys = 4M|r, in particular: \G'^ (4M|r)|/|GS(4M|r)| ~ 1.9. In the 
proton case, instead, as shown by the solid curve on the tight panel of fig. 11, is the S-wave that gives the main 
contribution, at the threshold: |GQ(4M^)|/|Gg(4M^)| ~ 0.2. 

Finally to have a comparison with data in the TL region, we consider the so called effective FF, G^s{(l^), cor¬ 
responding to the useful working hypothesis of a unique TLFF, that is: \G%{q^)\ = \G^{q^)\ = G^g{q^). Its 
expression in terms of the Sachs FFs follows by writing the e^e“ —> NN total cross section, obtained from 
eq. (33), as 


= <^pl(<7^) • [G2f(<7^)] 


47ra^ 

3(/2 



2A4 




(34) 


where (TpL(g^) represents the cross section in case of point-like fermions in the final state, which is obtained by 
putting G^ = G^j = 1 in the last expression of eq. (34). It follows that the effective FF is 


Geff(<7 


O'iVJv(g^) 

o-pL(g2) 



(35) 


Figure 12 shows the results for the proton (left panel) and the neutron (right panel) effective FFs together with 
all the available data. In case of the proton, left panel, the predictions, in particular the relativistic-uncorrected 
one, describe quite well the data in the high momentum transfer region, from 5^ ~ 7 GeV^ on, while they fail in 
reproducing the experimental G^g at lower q^, close to the physical threshold. Concerning the neutron effective 
FF, right panel of fig. 12, the predicted behavior does not agree with the available data that, however, cover only 
the near-threshold region. 




Figure 12. Effective FFs of the proton (left panel) and neutron (right) panel. The solid and dashed curves represent predictions 
with and without relativistic corrections. The vertical lines (red in the on-line version) indicate the physical threshold 
and the empty points are the world data sets form Ref. [8] and references therein. 


Finally, fig. 13 shows the modulus of the ratio electric to magnetic proton FF in comparison with the data. The 
isotropy-violation is manifest, having at the threshold a non-unitary value. The agreement with data, that favor 
a constant behavior at high q^, is quite poor, because, both results, corrected and uncorrected, have an increasing 
behavior, almost linear in q^. This is a consequence of the different high-g^ behaviors predicted for G^^{q^) and 
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both, in case of uncorrected results, where it is found oc ^ and G^{q^) oc \ see 

eqs. (C.17) and (C.19), and in case of the relativistic predictions, given in eqs. (14), where G^{q^) oc [constant] 
and G\,i{q^) oc 



Figure 13. Modulus of the ratio electric to magnetic proton FF, solid and dashed curves represent predictions with and without 
relativistic corrections. Two incompatible sets of data are shown: the circles are from the BaBar Collaboration [46] 
and the squares from the Lear Collaboration [47]. The vertical line (red in the on-line version) indicates the physical 
threshold. 


It is just such a failure in predicting the perturbative QCD power-law, see Sec. 2.3, that precludes the possibility 
of drawing any conclusion about the asymptotic regions. 

3.4. Isotropy at the physical threshold 

Following the treatment given in Sec. 1.3, isotropy at the production threshold manifests itself through the identity 

G%{m%) = GU^M%), (36) 

for proton, N = p, and neutron, N = n. Moreover, being the Sachs FFs (independent) linear combinations of the 
isospin components, i.e., G'g’)^^(g^) = G%^M(.q^) G^^Miq^)-, Ihe identity of eq. (36) is equivalent to 

Gf"" (4M^) = Gli^ (4M|r). (37) 

As already discussed in Sec. 2.2, the combination of such isospin components, that represent our primary outcomes, 
to obtain proton and neutron Sachs FFs, has to be performed with some care due to the different orders in 1/Nc 
expansion in which they are computed. In particular, from the definitions of eq. (11) and having that both, the 
mass M and the moment of inertia A are 0\Nc\, we get 

G% = 0 [V°] , Glt = 0 [iV°] , Gl = 0 [iV-i] , Gl, = 0 [iV,] . 


It follows that the more reliable test bed, as necessary condition for the isotropy hypothesis, is the isoscalar 
identity of eq. (37). In other words, the violation of such an identity would imply anisotropy. Figure 14 shows 
real and imaginary parts of the four isospin components of the electric and magnetic FFs in the TL region, across 
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Figure 14. 


Real and imaginary part of electric (dotted curve) and magnetic (solid curve) of the isoscalar (a,b) and isovector (c,d) 
FFs in the TL region, across the physical threshold 9phys — 4M^, indicated by the vertical red line. 


the physical threshold. In every instance, and hence also for the isoscalar FFs, figs. 14a and 14b, the identity is 
violated, i.e., the curves do not cross each other at the threshold, which is indicated by the vertical (red in the 
on-line version) line. 

Let us consider in more detail the constraints imposed by the isoscalar equation. The TL expression of 
is obtained by following the procednre, described in App. B, that has been used to compute the expression of 
G%{q^) given in eq. (B.6). In particular, considering the same symbols, it reads 


Gii{q^) = - ^ ^ Re {A, H [(-g - 3M^)zi] {zj + 1) + A, H [(? - 3M^)zj] {zj - 1)} 

+ + 1 ) - - 1 )] . 

At the physical threshold, q = 2Miv, the isoscalar magnetic FF is 

m 

= {R:^jH\{-2Mn - 2M^)zj] (1 -h Zj) - R.^H[{2 Mn - iM^)zj] (1 - z,)} 

N 

+ i^ E + c,) + - c*)] , 


while the electric one, from eq. (B.6), 


G|(4M^) = ^ Re {^R,^H[{-2Mn - 2,M^)zi] - R.^H[{2 Mn - 3M.)2,]} 


-f 


2 iTr 


1=1 
2 h 




Rc^e 


{2Mn-3M7t)cj 


S{yj) + Rle 


p* { 2 M ^— 


i=i 




It follows that the isotropy condition of eq. (37) becomes 


m 

'y ^ R,e H [(— 2tM]\[ — 3Mn^Zj] (1 + + 4A/ivA) — R^j [(2Afiv — (1 — 4AfivA)= 

h 

2 in Y,e{-Xj) [Ae^""”“®'^"’"^e(i/i)(l + Cj+ 4MjvA) -h - c* -f 4MjvA)] . (38) 

1=1 
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It can be interpreted as an implicit relation among poles (they appear in the argument of the H{z) functions, 
defined in App. A) and the corresponding residues of the function 6fit(r), that parametrizes the profile function 
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6(r), see eq. (17). 

It is a quite hard task to obtain the identity of eq. (38) from the beginning, i.e., as a condition which is automat¬ 
ically fulfilled by any parametrization. In fact, the possibility of using the definition of eq. (12) to relate directly 
the positions of the 6(r) poles to the properties of the chiral angle F{r), is prevented by the fact that such a 
relation holds only for real and positive values of r, while the poles Zj, j = 1 ,... ,m, lie in the r complex plane 
outside the positive real axis. In other words, by solving numerically the Euler-Lagrange equation of the Skyrme 
model, no information about the complex structure of the chiral angle F{r) can be accessed for r ^ (0, oo). 
Following the definition given in eq. (12), simple poles of 6(r) can be related to branch points of the chiral angle 
F(r). By considering eqs. (17) and (18), and assuming the coincidence between fit function and 6(r), we have 


h{r) 


foflt(r) 


F'{r) sin^ [E(r)] ^ y Rk 

2 tt^ r — Zk 


Such a differential equation for F{r) can be integrated and, by using the condition Rk = Oj it is 

fF{r) 


r' - Zk 


/ sin^ (F)dF = -27r^ '^Rk 

7f{0)—TT Jo 

2F{r) - sml2F{r)] = -8Tv^Y^RkZkl - + Zke-^^*''^>‘[Ei{3M„{zk - r)) - Ei{3M^Zk)] L (39) 


where Ei(z) is the multi-valued exponential integral function'*, that has branch points in 2 : = 0 and 2 = 00 . As a 
consequence, the function in the right-hand-side of eq. (39), besides the one at infinity, has m branch points in 
each r = Zk, with fe = 1, 2,..., m. A similar complex structure is expected for the chiral angle F{r), even though 
no explicit solution can be obtained due to the implicit nature of the left-hand-side expression. Figure 15 shows 
the analyticity domain of F{r) in the case where the branch cut of Ei( 2 ) is placed over the positive real axis. 
The cuts are obtained by adding to the negative real axis (negative because r appears in the argument of Ei( 2 ) 
with a minus sign) the points of the set {2fc}J!Li, that in the figure are organized in pairs of complex conjugates 
{cj, and real values with m = I + 2h, as in App. B. 

It is interesting to notice that, not only the behavior at the physical threshold, but the entire structure of TLFFs 
is intimately connected with the analytic extension of the chiral phase F{r) outside the positive real axis, which 
represents its natural domain. Moreover such an extension drastically changes the character of this function, 
because, by acquiring a non-vanishing imaginary part, it looses its ’’phase” nature. 


4. Conclusions 

A procedure to compute nucleon TLFFs, starting from integral representations of their SL counterparts, has been 
defined. Such a procedure consists in modeling the numerical solutions obtained for the nucleon electromagnetic 
currents in the framework of the Skyrme model of hadrons with particular analytic functions, whose Fourier 
transforms, not only, are well defined, but they also embody all the main theoretical features required, for the 
FFs, by first principles. 

The general form, for the analytic functions of the radius r, that describe the numerical solutions, is conceived to 


* The exponential integral function is defined as [44] 


/ OO 



t 


dt , 


the integration is in principal value for real and positive 2 . 
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have automatically the expected behaviors in the origin, r = 0, and in the limit r —>■ oo. 

The results for the nucleon FFs are analytic functions of or equivalently , which are real in the whole 
SL region and in the small portion of the TL region below the theoretical threshold while are 

complex elsewhere. Moreover, they also have the branch cut discontinuity oo), in the complex plane, 

as expected by assuming analyticity and unitarity. 

Once the analytic expressions for all nucleon FFs in the whole q^ complex plane are known, any quantity can be 
predicted without any further assumption or restriction. Indeed, the free parameters of the fitting functions can 
be fixed at any desired degree of precision, since the numerical solutions can be known with an arbitrarily high 
accuracy. 



Figure 15. Analyticity domain of the chiral angle F{r). The empty circles indicates the branch points and the lined bands are 
the corresponding branch cuts, that have been chosen with constant imaginary parts. Each branch point for F{r) 
corresponds to a simple pole for the profile function b{r), as a consequence of the definition given in eq. (12). 

It is important to stress that TLFFs are purely relativistic quantities, they can be defined only in the framework of 
the Quantum Electro-Dynamics (QED) and describe the vertex 'yNN, where a virtual photon produces a nucleon- 
antinucleon pair. In light of that, the only definition of an analytic continuation in TL q^ represents by itself the 
first relativistic extension of the starting FF expressions given in eq. (11). In spite of that, since such original FF 
expressions have been obtained for static nucleons, i.e., nucleons in their rest frame, a procedure has to be defined 
to extend FFs even at relativistic momenta. In the SL region we have adopted the approach described in Ref. [42] 
and the resulting, relativistically corrected FFs are shown in eq. (13). For TLFFs a modified methodology [20] 
has to be used to account for the behavior in the unphysical region, where FFs remain unchanged, and at the 
production threshold q^ — gp^ys, see eq. (14). Even though such a procedure does not reproduce the asymptotic 
behaviors expected from perturbative QCD, as it is also discussed in Ref. [42], the relativistic corrections improve 
the agreement with the dipole FF. 

The predictions have been compared with all the available data in SL and TL region. The fair agreement that is 
obtained, in most of the cases, in not negligible g^ intervals, appears as a quite encouraging achievement, since 
these results are based on a microscopic model which contains only pion fields and has no free parameters. 
Particular attention has been paid to the TL physical-threshold behavior, to verify if the non-trivial relationship, 
that exists between the predictions for the electric and magnetic FFs, reproduces the identity — 

expected in case of isotropy and non-singular Dirac and Pauli FFs. We observed that the complex 
equality is not fulfilled, i.e., the two independent equations for the real and imaginary parts are not verified. 
By considering the moduli, their differences are partially compensated, nevertheless the isotropy violation at the 
production threshold remains an important effect. 

In the TL region, especially nearby the threshold, the obtained values of the effective proton and neutron FFs 
are too small with respect to the data, while, at high g^, especially the relativistic-uncorrected ones, appear in 
better agreement with data. However, as already stated in Sec. 2, the Skyrme model represents an effective 
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approximation of QCD at low energy so that, pushing its predictability at high-g^ goes beyond the aim of model. 
Moreover, even the failure of relativistic corrections is expected, because it is well known that such corrections 
do not reproduce the perturbative QCD power law of FFs, which describes quite well the data. 

Since the present model contains only pions, the most natural improvement would be the inclusion of vector 
mesons p and uj as gauge bosons of a hidden symmetry [48]. This will entail additional degrees of freedom, i.e., 
further profile functions in terms of which parametrize the nucleon electromagnetic currents and hence the FFs. 
Finally, this procedure appears quite suitable to be applied to any other effective model of low-energy QCD, which 
allows to compute nucleon electromagnetic currents and then SLFFs as their Fourier transforms. 
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A. The integral representation of Ei{z) 

The function H{a/ 3 ) is defined through the integral representation 

rcxi —ar 

that, with a,P £ C, converges if Re(a) > 0 and P ^ ( 0 ,oo). Such an integral representation and hence the 
function H{aP) depend on a and P only through their product, indeed, by putting w = a(r+ /?) and z = a/ 3 , we 
have 


ntx -w 

H(z) = e’‘ / - dw = E\(z). 

Jz 


where Ei{z) is the exponential integral (Expintegral) function [ 44 ]. A series representation of Ei(z) can be 
obtained by integrating the well known expansion in the origin of its first derivative, i.e., 


dEi e-" 1 ^ (-2)'= 1 ^ 


(llJjl C i V "T Z) i V "T 

dz z z ^ k\ 

fe=0 


2 ‘ k\ 

fc=i 


Indeed, since the series converges uniformly, it can be integrated term by term as 

{-zf 


E,{z) = C-Hz)-Y, 


kk\ 


where C is the integration constant whose value is obtained by considering the limit 2 —>■ 0 , as 


C = lim 

z—yQ 


£:i(2) + ln(2) + ^^ 


lim 

z —^0 


j: 


- dw + ln(2) 

w 


7 is the Euler-Mascheroni constant [ 44 ]. In light of these results, the function H (2) has the representation 


H{z) = e" 


-7 - ln(2) - 


i-zr 


kk\ 


-7 - ln(2) + ^ 




kk\ 


It possesses the same properties of E\{z), i.e., it is analytic in the 2 complex plane with the cut (—oo, 0 ), for 
[ arg(2)[ < TT. It is real for 2 G ( 0 ,00), which is the interception of its analyticity domain and the real axis, so that 
it fulfills the Schwarz reflection principle 


H{z*) — H*{z), V2^(—oo, 0 ). 

B. The branch cut in the complex plane 

To study how the logarithmic cut of Ei{z), in the variable 2, evolves in the variable Q, see eq. ( 24 ), we consider 
the following cases 


2± = (±iQ — 3 MTr)zj , with: Zj = Xj + it/j , 
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where Xj and yj are the real and the imaginary part of the pole Zj. The variable Q ’’feels” the cut when z± crosses 
the negative real axis, i.e., when Im( 2 ;±) changes sign and Re( 2 ;±) < 0. The value Qo at which the imaginary part 
of z± vanishes and the corresponding real part, are given by 


Qo = ± 


SM^ryj 


Re( 2 :° ) = z\ 


+ y]) _ SM^lzjl"^ 

Xj Xj 


(B.l) 


where £ R stands for the value of z± corresponding to Q = Qo- 

From eq. (B.l) follows that: Re(z5-) < 0 when Xj > 0 hence, having Q > 0 by definition, there must be i/j > 0 
or yj < 0 depending on ±iQ. Moreover, if yj > 0 {yj < 0) the crossing is from above (below) the cut and so it 
requires the imaginary part to be increased by —2tt (+27r), having Ei{z) a — log(z) term, see eq. (23). 

Finally, all these considerations can be summarized in the compact expression 


Esl {z±) = El {z±) ± S{xj) 6{±yj ), (B.2) 

where the Heaviside 9 functions select the above conditions and the symbol Esl stands for an Expintegral 
corrected in case of SL momenta, i.e., Q € (0, oo). 

A similar study can be done also in the TL region. However, as already discussed, in such a region FF values, at 
a given above the threshold and (2M^)^ for isoscalar and isovector FFs respectively), are obtained as 

the limits 


G(g^) = Urn G(g^ + ic), 9 ^ > (3 - I)^M^ , 

e -»0 + 

where the symbol G stands for one the four nucleon Sachs FFs and / = 0,1 is the isospin. It follows that a generic 
value of 5 > (3 — I)M^, is understood as g + irj, with rj —>■ 0^. To obtain TLFFs from the expression of eq. (24), 
we have to make the substitution Q —> i{g + irj), and hence the arguments of the Expintegral functions become 

(±iQ — 3MT^)zj —>■ (=P 5 — 3M^ =p = (T<? ~ 3M-^)xj ± igyj — i [{SMt^ ± g)yj ± rjXj] . 

The imaginary part, being g and rj positive, vanishes only in the case of ’’lower sign”, at 

g = go = 3M^ - rj—. (B.3) 

Vj 

Since in case of G%, g > 3M^, the ratio Xj/yj must be negative, the pole Zj lies either in the second or in the 
fourth quarter of the 2 ; complex plane. The corresponding real part is 


U .|2 

Re[(qo - 3M^ + irf)Zj\ = (go - 3MT,)xj - rjyj = -g '-E- . 

Vj 

A correction has to be considered only if such a real part is negative, i.e., yj > 0 and, since Xj and yj have 
opposite sign, the only possibility for a pole to generate a correction in the TL region is that it must lie in the 
second quarter. When g —>■ g^ , of eq. (B.3), the imaginary part 

Im[(g - 3M^ + ig)zj] = -{3M„ - g)yj + gxj ~ (g - 3M„)yj 

vanishes as hm^_^g^+(g — 3M^)i/j = O'*" because yj > 0. So, following the previous argument, the imaginary part 
of the Expintegral will be increased by — 27r, hence we can define 

Etl [(g - 3M^)2 j] = El [(g - 3M,r)2j] - 2-Ki9{g - 3Mn)9{-Xj)9{yj). (B.4) 
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These corrections are crucial because, as we will see in more detail, they generate the desired complex structure 
for the FFs. 

Having real polynomials with only simple zeros, the poles of 6fit(r), {zj}JLi, can come either as single real negative 
values or in pairs of complex conjugates {cj, c*}j=i, and hence U {cj,Cj}j^i, with: 

I + 2h = m. Moreover, from the definitions given in eqs. (19) and (21), the residues have the same properties of 
the corresponding poles, i.e., 


Aj € K , j - 1 , 2 ..., 1 , . 

3 

In light of this, using the function H{al3) of eq. (23), in particular its property: H{z*) = H*{z) (Schwarz reflection 
principle) and including the branch cut corrections of eq. (B.2), the expression of G%{Q) in the SL region, given 
in eq. (24), can be simplified as 


4 ' 4 

G|(Q) = [(*'3 “ ^ I] Im [Rc^H[{iQ - 2.MYcj\ - R^^H[{iQ - 3M,,)’c,]} 

^ t=i 




+ — ^^(xjORe 
i=i 


e[Q-‘iM^^]e{y,)Rc,e 


^(iQ-3M^)cj 


j +0(Q + 3M^^)e{-yYRo,e 




The TL expression of eq. (24), accounting for the corrections of eq. (B.4), becomes 


G%{iq) = - — ^ Re {[(-g - - Ro^H[{q - iM^Cj]} 

1=1 

9 ‘ 

— -Y. - SM^rj] -H[{q- SM^Vj] } 

1=1 

+ ^e(g - 3MYe{-xY • (B.6) 

1=1 

While the SL G%{Q) is real, G%{iq), containing complete H functions (not only their real or imaginary parts), 
could have a non-zero imaginary part. 


C. The asymptotic behavior 

The integral representations (Fourier transforms) of FFs given in eqs. (11) can be classified into two species 
depending on the order of the spherical Bessel function, i.e.. 


go{Q) = 

coo 

/ f{r)joiQr)dr, 

JO 

(C.7) 

si(Q) = 

coo 

/ f{r)ji{Qr)dr, 

Jo 

(C.8) 


with 


3o{x) 


sin(a:) 

X 


MQ) 


sin(a:) cos( 2 ;) 

x"^ X 


djojx) 

dx 


and /(r) represents the profile function, which is regular in the origin and vanishes exponentially as r —l oo, in 
particular 


f{r) 


oc r* , 

r—>-0 


/(r) oc 

r^oo 



(C.9) 
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with /, /i e N and fj, > 0. The behavior in r = 0 is crucial because it determines the asymptotic trend, as Q —>■ oo, 
of the functions gi,o{Q)- The profile /(r), which is known only numerically, is parametrized as 


fir) 


^m(r) —fir 

B„{r) ^ 


where Amir) and -B„(r) are the real polynomials 


m n 

Amir) = '^akr'", Bn{r) = '^bkr’° , 

k=l k =0 

with n — m = h, and a; 7 ^ 0, am yf 0, &o = 1, 6 n yf 0, in order to follow the behaviors given in eq. (C.9). 
Assuming the zeros {zk}k=i of -Bn(r) to be all simple, with Re( 2 fc) < 0 and also Amizk) 7 ^ 0, Vfc G {1, 2,..., n}, 
the polynomial part of /(r) can be written in terms of the Mittag-Leffler representation 


fir) 


fc=l 


Rk 

r - Zk ’ 


Rk = Res 


Amir) 

B„{r) 


,r = Zk 


where Rk is the residue of the simple pole Zk- The two polynomials have n + m — Z + 1 real degrees of freedom, 
while the Mittag-Leffler representation has 2n>n-|-m — Z-|-l free parameters, namely the complex zeros Zk and 
residues Rk, with A: = 1, 2,..., n. Indeed, having Bn{r) real coefficients, both, zeros and corresponding residues 
of the ratio, can be either real or pairs of complex conjugate, hence the real degrees of freedom are only 2n. 
Summing up the Mittag-Leffler series we get 


n ^ — 1 -j' 

Rk _ r^jr 

^^r-Zk Bn{r) 

where the coefficients aj contain zeros and residues. To have, at numerator, a polynomial of m < n degree with 
a zero of order I in the origin we should impose 

aj = 0 , V j G {0,1,..., Z — 1, m -I- 1, m -I- 2,..., n — 1} , 

lvalues n — m —lvalues 


these are Z-l-n —m — 1 constraints. It follows that the number of degrees of freedom reduces to 2n — (Z-l-n —m — 1) = 
n -I- m — Z -I- 1, which coincides with that of the ratio of polynomials. In particular, the zero of order Z in r = 0, 
first condition of eq. (C.9), implies 


k=i * 






/‘'■'Ho) = 0 


C.^±^=0, 

fe = l 




(C.IO) 


fc=l 
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In general, the Cj constants are real, being sum of real and/or pairs of complex conjugate numbers. 
The integral representation of eq. (C.8) can be put in a form similar to that of eq. (C.7) as 


^OO 

9i{Q)= f{r)ji{Qr)dr 
Jo 


9i{Q)=- 


jojQr) 

Q 




OO ^ 

0 Q 


/’OO 

/ jo{Qr)f'{r)dr 

Jo 


1 

S'i(Q)=gy^ jo{Qr)f{r)dr. 
Assuming Z > 1, we define g{r) = /(r)/r, with 


(C.ll) 




and the Mittag-Leffler representation 


n-|-l = n' _ 

9{r) = e-^^ E — 


hence the derivative 


/'(r)=rsr'(r)+fl(r) = e E 

C.l. The function goiQ) 


gr — 1 
r - Zk 


(C.12) 


The analytic expression of go{Q) can be obtained by integrating the representation of eq. (C.7). In particular we 
have 




sin(Qr)e-'^’~ A J^dr = - 

r E—/ 


r - Zk Q 


1 Rk r 

^ fe=i 


sm{Qr)e 


r - Zk 


~dr — 


f 

Jo 


sin(Qr)e 


'dr 


where the first integral can be computed in terms of the Expintegral function, while the second, which does not 
depend on Zk, is equal to the arctangent of Q/g. Hence we have 


fc=l 


’e-(^-'«)"'=Ei[-(M - iQ)zk] - + iQ)zk] 


2i 

\t^J _ 


The last term vanishes with I > 1, i.e. if /(O) = 0. In general, the profiles that appear in the integrals of eqs. (11) 
are bounded in the origin and hence, thanks to the factor of the differential d^f, the function /(r) has always 
a zero in r = 0 of order I > 2. It follows that the expression of go{r) becomes 


k = l 


e-{'^-'«)""£;i[-(p - iQ)zk] - + iQ)zk] 


2i 


The asymptotic behavior of the Expintegral function can be derived from the expansion [45] 

-2 N 


i=o 
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We consider go(Q) in the limit Q —>■ oo 


k=l 3=0 


1 oo ^ p 

^ j=0 fe=l 


1 


[-{^1-iQ)zhY+^ [-(/r + iQ)zfe]'’+^ 
1 1 " “ 


+ iQ)t+i - iQ)t+i 


_ ^ ■I/-' Irn [(m “ 

^ n ^ J-<^j+2 ^^ 2 , cr2y+i 




(m" + Q'^Y 


where the last snm starts from j = I — 1 because the coefficients Cj, defined in eq. (C.IO), are vanishing for j < 1. 
The imaginary part can be written in powers of Q so that 


5o(Q) = 

j=i-i 


YCj+2 

(^2 + Q 2 ),+r 


Int[i/2] 

E 


s = 0 


2s+ 1 

i + 1 


s+1^2s j — 2s 




(C.13) 


The highest power of Q in the numerator coincides with the maximum even number less or equal to j. Hence, 
for two even-odd consecutive values of j, the highest power of Q at numerator remains the same, while that 
at denominator increases linearly with j. This means that the terms at higher orders in j are higher order 
infinitesimals as Q —>■ oo. In particular, the j'-th term behaves as 


O 


g2Int(j72)-20'-|-l)j 


' O = O [qY 

O = O [Q-"] 

O = O [Q-®] 


j = 1 
i = 3 ’ 


(C.14) 


the dominant asymptotic behavior is given by first two terms, with j = 1 and j = 2. 
In the TL region, i.e. Q = iq, with g > 0, the expression of eq. (C.13) becomes 


9o(ig) 


OO 


= -E 


YCj+2 


Int[j72] 

E 


s = 0 


2s -f 1 

J + 1 


2s j — 2s 

q M 


(C.15) 


hence the asymptotic behavior follows the same power law of eq. (C.14). 

More in detail, once the order I (see eq. (C.9)) of the zero, that the profile function posses in r = 0, is known, 
also the asymptotic behavior in both, SL and TL regions, is obtained as 


{l-l)\a+r .,(,+i)/2W-i (^-l)!Ci+i(-l)<^+^)/^ 

(^2 + g2)! ~ Ql+l 


9 o{Q) < 


n(-i) 


1/2 


(m'+Q 2 ); + 1 


(/rCj+i — C1+2) + 


3 „ „,_ 2 i - C 1 + 2 ) 


Q '+2 




go{iq) —^ < 

q—¥oo 


{g? - q^y 
l\ 

{9-^ - q^y+^ 


Qf^ + l 


I odd 


I even 


I odd 


(C.16) 


1^—q^(/rC'i+i — C'i+2) + fYCi+iq‘ 


3^ 1-2] l\{fiCi+i — C1+2) 


ntZ +2 


I even 


The electric, isoscalar and isovector, FFs, eq. (11), are obtained through integral representations of type (C.7) 
with I = 2 and I = 4 respectively, and hence 


Ge{z) ~ , 

2 —>00 



Z —>-oo 


(C.17) 


where, the SL and TL limits are considered with z = Q and 2 : = iq, respectively. 
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C.2. The function gi{Q) 

The asymptotic behavior of gi{Q) can be achieved by the integral representation of eq. (C.ll) and the expression 
of the /(r) derivative given in eq. (C.12) as 


9i{Q) 


^ roo 

/ sin(gr)e“ 

^ fe=l •'0 


N' 


(/r - 


+ (l/zfe - g) 




-1 1/^ - 9 _ J_\ 

(r - Zk)^ r - Zk Zkr ) 

Ei[-{g - iQ)zk\ - (m + + iQ)zk] 

2i 

— iQ)zk] — + iQ)zk] arctan(Q/^) 

2i Zk 


The terms proportional to jj, cancel and that proportional to the arctangent, assuming Z > 1, is vanishing and, 
using z = —(fi — iQ), we have 


n' 




e^^^Eiizzk) + e" ^^ Ei{z*Zk) ^ e^i{zzk) - e" 


'^>‘Ei{z*Zk) 


2Q 


2iQ^Zk 


By taking advantage from the asymptotic series of the Expintegral function, the first term can be written as 


e^^>^E^{zZk) + e^'^'^E^{z*Zk) 1 , 


2Q 


1 1 
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-E 
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E 
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while for the second term we have 

e^^kEYzzk)-e^’‘^kE,(z*Zk) 

2iQ^Zk 

The complete expression is then 


-E 


Int[3/2] , I 1 

s! V r I 2t T 1 

^ ' s + 1 


23+2| | 2 s +2 

=0 \^\ t=0 


(- 1 ) Q g 


gYQ)- E 


s!C' 


3 + 1 


(/i2+Q2)3+1 

s!C '+2 


Int[{s + l)/2] 

E 

t=o 

Int[s/2] 

E 


2t 

5+1 

2i + 1 
5+1 


(-1) Q 
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(-1) Q g 


where the constants Ct are defined as those of eq. (C.IO), but for the residues Dk, i.e. Ct = 'YYYk=\ Ekjz\ and the 
lower limits of indexes s account for the behavior of /(r) at r = 0. The TL asymptotic behavior can be obtained 
from the previous expression, by setting Q = iq, with q —>■ oo, i.e. 


g\ (*<?) 
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s=Z —1 
oo 

+*E 
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slC's^l 
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The leading terms are 


ffi(Q), 


Q'+i 

(_l)(i+i)/2(^ _ i)!(; + 1 _ 5 i.0(qVi - mC;) 


rr^+1 




a-i)!(^ + i-5i,o(Q'+i-Mc;) 




I even 


I odd 


I even 


I odd 


(C.18) 


The magnetic, isoscalar and isovector, FFs, eq. (11), are obtained through integral representations of type (C.8) 
with an additional factor Q~^, or {iq)~^, and Z = 3 in both cases, it follows that 


z—¥oo 


Gm{z) ~ 2 

z—foo 


(C.19) 


where, as in the go case, SL and TL limits are considered by setting z = Q and 2 = iq, respectively. 
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